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1 Introduction 



In this paper we will analyse an SU{2) Yang-Mills field configuration which is periodic in 
3-space and tends to a pure gauge in both t = ±00. This configuration is a solution of 
the classical euclidean equations of motion and has infinite action and topological charge 
( finite within one unit cell). The configuration arises naturally when one studies gauge 
fields on the spatial torus with twisted boundary conditions (t.b.c.) [1]. In particular, we 
fix the torus to have equal period in all three directions and a twist vector m = (1, 1, 1). 
The presence of the torus breaks down the 5*0(3) rotational symmetry down to the cubic 
group. If we fix, without loss of generality, the length of the spatial torus to / = 1, we 
may write the boundary conditions 

A^(a; + ii) = criA^{x)ai (1.1) 

where cij are the Pauli matrices and Cj the unit vector in the i— direction. The choice of 
the Pauli matrices as the twist matrices can be regarded as a partial gauge fixing of the 
problem. The remaining group of gauge transformations must satisfy 

fl{x + ii) = ain{x)ai (1.2) 

This is not, however, the most general internal symmetry transformation of our 
problem. One can perform a transformation 

A^(x) -^ A'^{x) = n^^'\x)A^{x)n^'^^{x) + in^'\x)d^n^'^^{x) (1.3) 

with 

n^^'\x + e,) = z, aM'\x)a, (1.4) 

and Zi = ±1. This transformation preserves the b.c. Eq.(l.l) and leaves the Yang-Mills 
action invariant. Nevertheless, ii Zi ^ 1 for some i = 1,2,3 the transformation is not, 
strictly speaking, a gauge transformation since it modifies the values of the Wilson loops 
associated with non-contractible loops (Polyakov loops): 

W{^) -^ W'{^) = zT^^'^ W{-i) (1.5) 

where a;, (7) is the winding number of the loop in the i— direction. We will call these 
transformations singular gauge transformations and the group of these transformations 
(modulo ordinary gauge transformations) is isomorphic to Z|. As we will see later, these 
transformations act non-trivially on our gauge field configuration and thus we have 8 
configurations rather than one. 

The configuration with smallest action (S=0) with our boundary conditions is of the 
pure gauge form: 

A^{x) = i n^^'\x)d^n^'^\x) (1.6) 

with fl^^\x) satisfying Eq.(1.4). However, it turns out [2] that there are only two gauge 
non equivalent configurations of this type, those with Zi = 1 for all i, and Zi = —1 for 
all i. In the Aq = gauge these configurations are constant in time and equal to the 
minima of the potential energy J2iTrBf. The problem becomes that of a Z2 symmetric 



potential with two non-symmetric minima. The order parameter distinguising these two 
is the Polyakov loop winding the torus once along each direction {uji = 1). Just as in the 
case of the X(f)* potential (and negative mass squared), we can investigate the solution 
which interpolates between the two minima as time goes from — cxd to +00. This is an 
instanton solution and gives rise to the leading weak coupling contribution to the energy 
splitting between the Z2 symmetric and antisymmetric states. 

These instanton configurations are precisely the solutions which we are after. They 
have half integer topological charge and are precisely the minimum action configurations 
in a 4-torus with infinite length and with t.b.c. both in time and space. The twist in time 
rioi = 1 eliminates the constant minima solutions. 

The existence of these configurations was proved by Sedlacek [3] . In a previous paper 
[4] we investigated the behaviour of these solutions by using the lattice approximation 
and the cooling method [5]. There it was shown that the total action and energy of the 
lattice minimum action configurations scale towards a continuum value. The value of the 
continuum action is very close to 47r^, the absolute minimum for a topological charge Q of 
1/2. Indeed, this minimum is only attained by a self-dual or anti-self-dual configuration, 
so our configurations must be very approximately, if not exactly, self-dual. If we perform 
a parity transformation to this self-dual instanton configuration we get an anti-self-dual 
anti-instanton configuration with Q = —1/2. 

The purpose of this paper is to report a more extensive and accurate analysis of the 
lattice minimum configurations in order to answer some questions which were not settled 
by our previous paper. Apart from giving additional support to the question of scaling to- 
wards the continuum and self-duality, we have shown that there are indeed eight different 
gauge-inequivalent configurations (four instantons and four anti-instantons) which can be 
obtained by acting with the group of singular gauge transformations on one of them. All 
these configurations are non abelian and cubic-symmetric. We have obtained functional 
expressions for the various physical quantities, including the vector potentials and field 
strengths in a suitable gauge, which describe their qualitative features. Satisfactory quan- 
titative description of the data has been obtained by using these functional expressions 
and the first few terms in a Fourier expansion of the functions involved. We also report 
the steps that we have taken in the direction of finding an analytic expression for the 
solution. Although our attempts in this respect have not achieved the ultimate goal, we 
have explored several paths involving some lenghthy and non-trivial manipulations, and 
our results can be of great help for future investigations. 

Our solution, being non-abelian, is very different to the abelian solutions [6] which 
are known to exist for other twists, other values of the topological charge and additional 
periodicity in time. However, it is much more similar to the instanton solution on S3 x K. 
In a recent paper [7] this solution has been used to investigate the ^-dependence beyond 
the steepest descent. It is in this spirit that our solution is physically relevant. The reader 
is addresed to Ref.7 and the earlier Ref.8 to see the role of our instanton solution for the 
Yang-Mills dynamics, in the weak coupling limit. 



2 Analysis of the data 

Our strategy has been explained in Ref.4. We consider an A^^ x A^^ lattice with twisted 
boundary conditions. We use cooling to obtain configurations with smaller and smaller 
action. We stop once the value of the action has changed by less than 10~^ over the last 100 
cooling sweeps. To perform the analysis reported here we have used the configurations 
which we employed in Ref.4 together with new data which extend over larger lattices 
{Ns,Nt) = (11,29) and (15,29). We have little to add to the results of Ref.4 concerning 
integrated quantities. The new data just give extra support to the conclusions presented 
there. The value of the action S is 39.234 and 39.347, for (11, 29) and (15, 29) in agreement 
with our extrapolation formula S = 47r^ — 29. 7a^ (with a = l/Ng). The degree of self 
duality is now X = 2 g^~ ^ = 0.0045 and 0.0033 for the same two lattices. We notice 
that the errors introduced by the lattice approximation are smaller than 1% for integral 
quantities. This can be explained by the fact that they are 0{a'^). In the case of the 
vector potentials and the colour electric and magnetic fields at a given point A1{x), E'^{x) 
and Bf{x), we expect errors of order a: typically then 10% for our values of A^^. One can 
estimate more precisely the size of these errors by several methods. For example comparing 
E^ and Bf at the same point the difference is never greater than 0.14 ( representing 3% 
of the magnitude at those points). On the average a difference of order 0.04 makes the x^ 
per degree of freedom of the comparison of order 1. Other estimates of the errors can be 
deduced by comparing values for different Ns or by using different lattice approximants 
to the continuum value. The resulting errors depend on the quantity under consideration 
and on the value of t, but stay within a factor of two of 0.04 in all cases. We have chosen, 
thus, the value of 0.04 as a typical error which we will be using in all x^ fits of the data. 

2.1 Gauge invariant quantities 

We start analysing the structure of our solution by looking at gauge invariant quantities. 
Consider first the colour electric and magnetic fields E"-, B^. They can be regarded as six 
vectors in the 3-dimensional colour space. Gauge transformations amount to rotations in 
this space. Thus, the gauge invariant quantities are the scalar products and moduli of 
these vectors Mij = Ei-Ej. Due to the approximate self-duality the colour magnetic fields 
give the same information within errors. Before stating the results of our analysis let us 
mention one technical point. Our interest is centered on the values of these quantities in 
the continuum limit. To extract these values one can use various lattice quantities all of 
which differ by an amount of order a. Nevertheless, in checking certain symmetries, it is 
essential to choose lattice observables which respect the lattice symmetry. For that reason 
we have extracted the colour field strength at each lattice point, F^j^(na), by averaging 
the four plaquettes attached to that point in each /xz/ plane: 

F'^,{na) = -^Tr{-ia, [P^,(n) + P._^(n) + P_.^(n) + P-^_.(n)]} (2.1) 



i/—fi 



n 



where, for instance, P,^_^(n) stands for the following product of lattice links: P,, 
XJiy{n)\JJ^{n — fi + u)\j\{n — /i)U^(n — ji). The previous quantity transforms in the right 
way under gauge transformations and is symmetric under cubic transformations around 
the point na. 



The main features of the solution are the following: 

a) At each time value there is a point of maximum energy density, which we will call 
the spatial center of the instanton, whose coordinates are time-independent and we 
will choose them as a; = 0. Furthermore there is a time where the energy is maximal 
and we will choose it as the origin of time, t = 0. Our solution is cubic invariant 
with respect to the spatial center and invariant under time reversal. The spatial 
center is not one of the lattice points, but it can be located by interpolation. The 
same is true for the origin of time. 

b) For all times and at the spatial center of the instanton one has Mjj ex 5ij. This 
orthogonality implies that the solution is not abelian. Away from the center the 
vectors E^ cease to be mutually orthogonal but they are never collinear. 

c) We have fitted the moduli and scalar products at t = to the first terms in a Fourier 
expansion. Our results are 



Mil = cos2(7ra;) {15.00 + 20.02 [0(y) + 0(2)] 

-3.12 cos2(7r|/)cos2(7r2) + 0.128 sin'^{2Tiy)sin'^{2Tiz) 

-0.18 [cos^{Txy)sin^{2'nz) + cos'^{TXz)sin^{2ny)\] (2.2) 

+32.74 0(i/)0(^) 

-sin^{2Tix) { 0.232 + 2.25 [cos'^{'ny) + cos'^{ttz)\ 

+0.80 cos^^Tiyyos^^-Kz)} 

M12 = -1.608 [1-0.433 cos(27r2)] X(a:) x(t/) (2.3) 

where 

(t){x) = cos'^{nx) - 0.082 sin'^ {27rx) 

x{x) = sin{2iTx) + 0.086 sm(47ra;) 

The quality of the fit is very good ( see fig.l for example). The remaining components 
are just obtained by the apropiate replacement of the indices and variables in the 
previous expression, in a way which is consistent with cubic invariance and parity 
invariance. Notice that Ef vanishes when Xi = ±1/2 and Xj = ±1/2. Notice also 
that the solution seems perfectly smooth. 

We turn now our attention to non local gauge invariant quantities: the Polyakov loops. 
Due to the boundary conditions these variables are defined as 

n^{x) = Texp{i f A^dx"} r^(x') Texp{i f A^dx"} (2.4) 

where 7^(a, b) is a straight line in the positive /i direction starting at a and ending at b, 
x' is the border of the torus patch, Texp is the ordered exponential and r^(a;) the twist 
matrices. On the lattice these quantities are simply given by the ordered product of the 



/i-links corresponding to the path '~f^{x). These quantities transform hke Ej under gauge 
transformations and therefore the gauge invariant quantities are just the scalar products in 
colour space |Tr(0^(x)f2j^(x)) = X^^{x) and the traces themselves |Tr(f2^(a;)) = X^ 



X] 



Since our instanton evolves in time from one pure gauge Eq.(1.6) {zi = 1) to another 
gauge inequivalent one {zi = —1), and a representative of each class is given by fli{x) = 
±icri {i = 1,2,3), we expect Xij = —Sij, Xj = for large |t|. Indeed our results agree 
with this situation. As a matter of fact X^ can be fitted with a x^ = 0.005 per degree of 
freedom ( with an absolute error of 0.04) to the formula 

X^=l[m, (2.5) 

where rrii = m{xi) for i = 1, 2, 3 and ttiq = fno{t)- This is a very remarkable factorization 
property. The function m is parametrised as 

m{x) = cos{ttx){1 + Csin'^{7rx)) (2-6) 

where C = —0.196 and mo(t) is given in Fig. 2. Notice as well that the m functions are 
antiperiodic. Thus, by translating our solution by one period in any direction we get 
a new solution. In this fashion one generates 4 solutions out of one. Thus, we are not 
dealing with a unique instanton solution but rather with a family of solutions related by 
singular gauge transformations. 

2.2 Gauge dependent quantities 

In order to extract all the information from our numerical results we have to obtain also 
gauge dependent information. In particular, one would like to extract the gauge potentials 
themselves in a suitable gauge. This turns out to be feasible and our results are presented 
in what follows. 

First of all one must introduce a simple gauge fixing, which respects most of the 
symmetry of the problem. It is natural to select Aq = 0, a condition which preserves 
all spatial symmetries. To completely specify the gauge one has to fix the remaining 
time-independent gauge transformations. This can be done by fixing A{x,t = — oo) = 
or equivalently fli{x,t = — oo) = iai. This fixes the gauge completely including global 
gauge rotations. At t = +oo the gauge fields must coincide with a pure gauge Eq.(1.6) 
and in the Aq = gauge this gauge transformation is indeed given by the temporal 
Polyakov loop Q,o{x,t = — oo) . Notice that fl^ transforms precisely like Eq.(1.4) for 
Zi = —1. The exact form of flo{x,t = — oo) can be obtained in terms of the gauge 
invariant quantities |Tr(Oo(a^, t = — oo)) and |Tr(Oo(a;, t = — oo) ■ Cli{x,t = — oo)). 
These have been obtained from our data and the result can be summarised by giving a 
functional form which describes these data with a x^ = 0.01 per degree of freedom 

Q,o{x,t = —oo) = ■m{x)m{y)m{z) I + i 1 — Y\_''^'i — ^~ (2-7) 

V i \j\ 

where fi = f{xi) and f{x) = sin{TTx){l + Bcos^{'n-x)), with B = —0.172. Notice that f{x) 
is very similar to m{x + 1/2). Changing the value of 5 to C does not change significantly 
the quality of the fit. 



Now we turn our interest towards computing Ai{x,t) in the previously mentioned 
gauge. To extract these quantities from our lattice results implies some lengthy procedure. 
Apart from the typical order a errors we have to deal with some other sources of errors. 
Some errors arise through the gauge fixing procedure due to the finite value of Nt/Ng. In 
this case Aj is never a pure gauge and cannot be set to zero. Nevertheless these errors are 
for most purposes quite small and can be monitored by varying Nt/Ng. Our procedure to 
evaluate the gauge fixing on the lattice configuration is as follows. To evaluate the gauge 
fixed value of a lattice link U^(n), we choose a lattice path going through this link and 
starting and ending at some point P. This point is situated on the surface of minimum 
energy ( smallest time rit = —{Nt — l)/2) and the path is made of three pieces 71, 72 and 
73. The first and last paths are situated on the surface of time coordinate —{Nf — l)/2. 
7i joins P with P' = (— ^ *^~ ,n,), the point with equal spatial coordinates as the origin 
of the link U^(n). 73 joins P" = (— ^ *~ , n + jl) with P. 72 moves from P' along the 
positive time direction, then goes through the link in question, and then back to P" along 
the negative time direction. We complete our gauge fixing by performing a global gauge 
transformation ^{P) which rotates the spatial Polyakov loops i7j starting at point P to 
the values iai. Summarising these transformations we write down the expressions for the 
gauge transformed link U' (n): 

u;(n) = niP)v,,v,,v,,n\p) (2.8) 

Indeed, if the spatial links at the smallest time are a pure gauge configuration, U-^^ and 
U^3 can be gauged to /. And since in the temporal gauge U^j = U^(^)) we see that 
U|^(n) is really our gauge fixed link. 

In addition, the position of P and the choice of 71 and 73 are irrelevant. However, as 
mentioned previously we do not exactly have a pure spatial gauge for —{Nt — l)/2 and 
the actual choice of 71 and 73 is relevant. To preserve the fact that U' (n) (Eq.(2.8)) is 
a gauge transformed of U^(n), one must choose P, 71 and 73 in a predetermined way for 
all lattice links. Our choice is as follows. The path 71 is obtained by running first in the 
1 direction, next in the 2 direction and finally in the 3 direction. We finally obtain: 



V'Jn) = fl{P)L{n)\J^{np{n + fi)n\P) (2.9) 



with 



«2-l 



L(n) = n Ui(P + giei) ■ n U2(P + n[e, + ^262) • (2.10) 

<ji=0 52=0 

"3-1 "o^-l 

Yl U3(P + n[ei + 71262 + qsis) ■ J| Uo(P + n[ei + 71262 + n'^h + goCo) 

<?3=0 <jo=0 

where n' = n^ — P^ and P^ are the coordinates of point P. 

Once the data are gauge fixed we can extract the lattice approximants to the continuum 
fields. The potentials A'} are taken to lie on the mid-point of the corresponding links. The 
magnetic fields Bf, are computed by the plaquette averages mentioned before. Our results 
show the same qualitative features for both vector potentials and magnetic fields. These 
features are described by the following functional form which fits the data: 
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Ol = qi{t) cos{ny) cos(7r2;) (I + q2{t)cos{2TTx)) 
(l + qsit) cos{2Try)) {I + qs{t)cos{2Trz)) 

Of = (1 + ai(t) cos{2nx)) (1 + a2it)cos{2nz)) 

[a^it) sin{'Kz) cos{'kx) {1 + ai{t)cos{'Ky)) (2-11) 

+a^{t) cosijxz) sinijix) sin{2T[y)] 

Of = (1 + ai{t) cos{2ttx)) (1 + as (t) cos (27ri/)) 

[—a^it) sinijiy) cosijix) (1 + a4(t)cos(7rz)) 
+ a^{t) cosijiy) sinijix) sin(2Tiz)] 

where O" stands for either 5", E^ or A", but with different values of the parameters qi 
and Oj. The remaining components of O" can be obtained by cychc permutations of the 
indices 1,2,3 in Eq.(2.11). The values of the parameters appearing in the expression are 
given in tables 1 to 4. 

The previous functional form describes satisfactorily the data with a maximum x^ per 
degree of freedom equal to 4 for A" and 6 for B'}. It can be used to study the main features 
of the solution. Another interesting aspect of the formula is to serve as a guide for the 
construction of an ansatz which could lead to an analytic expression of the solution. In 
the following section we will describe the steps that we have taken in this direction. 

3 Ansatz 

In Eq.(2.11) we have given an analytic expression which describes our numerical results 
within errors. It is clear that the solution is not as simple as that expression, but we want 
to explore the consequence of the fact that it shares with it the same behaviour under the 
symmetry operations of the system. Thus, we will assume the following general form for 
the solution (for all t): 

At = 6tQ^'\t,x„Xj,Xk) + {e'''''fQ^^\t,x„Xa,xt) + e'^'Q^^\t,x„Xa,xt) (3.1) 

where Q^^\t, x, y, z) = Q''^~^\t, x, y, z) is symmetric under the exchange oiy and z, even in 
all the three variables, antiperiodic in y and z, periodic in x. Q^'^\t, x, y, z) is odd in x,|/, 
even in z, antiperiodic in x^z and periodic in y. Q^^\t, x, y, z) is odd in the third variable 
and even in the other two, antiperiodic in x,z and periodic in y. The previous expression 
implies that the solution is cubic symmetric with A1{x) and B^{x) transforming in the 
Ti Ti representation of this group (Ti is just the spin 1 representation of rotations). 
In addition, the behaviour under translations by one period is just what is required by 
the twisted boundary conditions. The previous properties are more restrictive than what 
can be deduced by these symmetry properties alone. More precisely Q^^\t,x,y, z) could 
contain an additional term Q'^^^\t, x, y, z) which is odd in all three variables, antiperiodic 



and antisymmetric in the two last variables y,z. To test the presence or absence of this 
term we have added to A\ in the parametrization Eq.(2.11) a term 

g4 sin(2Tix)sin{ny)sin{nz) {cos{2ny) — cos{2nz)) (3.2) 

and refitted our data for A and B. The presence of a such a term improves the value of 
the x^ by a factor 2/5, but the size of q^ and thus its contribution is of the order of the 
estimated size of the 0{a) errors. Thus, we will keep an open attitude and explore both 
possibilities. 

Another additional information which might be crucial for finding the analytic expres- 
sion is the behaviour under P ■ T. Both parity and time reversal map instantons into 
anti-instantons, but one could ask whether the product of these two leaves our solution 
invariant or not. The implications of this invariance are complicated by the future-past 
asymmetry of our gauge choice. Invariance under P ■ T implies 

- A,(-f, -t) = nlAiix, t)Oo + t^ldi^o (3.3) 

where fio is the gauge matrix which describes the pure gauge at t = oo. The matrix ilo 
was parametrised by Eq.(2.7) and indeed the form of the axis of rotation //|/| is precisely 
consistent with the restricted form where Q*-^-* = Q^^~^^ even in all three variables. The 
necessary and sufficient condition for the gauge field at oo to be of the restricted form 
q(i-) = is that fi is a function of Xi alone. 

Expression (3.3) confiicts nevertheless with the requirement Q^^^^ = since in principle 
a rotation by fio does not preserve this constraint. If we impose that a general rotation 
around //|/| should preserve Q^^~^ = we must have 

Q^'\t,x,y,z) = fiy)Q^'\t,x,y,z) (3.4) 

Q^'\t,x,y,z) = f{z)Q^'\t,x,y,z) (3.5) 

where Q*-^-* (t, x, y, z) and Q*-^-* (t, x, |/, z) are symmetric functions of the last two arguments 
y,z. We have verified that these additional constraints are reasonably well satisfied by 
our data. Thus, we arrive at a restricted parametrization consistent with our data: 

An^, t) = Pi 7^ + P2 \fm - ^) + P3 ^auh (3.6) 

I/I l/r 

where p|, pg ^"^^ P\ must be even functions of the three variables x^y^z. p\ and p^ ^-re 
periodic in Xi and antiperiodic in Xj ^ Xi, and pg is antiperiodic in Xi and periodic in 

The fact that this form is not a mere consequence of the symmetry properties im- 
plies that, when requiring self-duality of such a solution, new equations will arise which 
guarantee the form to be preserved. 

To start with, let us perform a time independent gauge transformation and write 

A[{x, t) = vt^^'^Ai{x, t)nl'^ + in^^'^d,vtli/2 (3.7) 

This new vector potential complies also to the form Eq.(3.6) in what respects to the 
behaviour of the p' functions under parity. In addition it must satisfy — A^(— x, — t) = 



A.'-(x, t). This implies that p^ and P2 are odd in t and p^ is even in t. We have exphcitly 
checked these properties in our nummerical data and found agreement within errors. In 
fact, given self-duahty it is enough that p/ = P2 = at t = to guarantee the appropiate 
behaviour for t 7^ 0. At t = we have been able to fit both Aj and Bj with p^ alone. 

Imposing now self-duality we obtain six equations rather than three. The first set of 
equations follows from the requirement that the p^ should have the required behaviour 
under parity. The equations are: 

e^^\V,p:' + T^'^p','p:')=0 (3.8) 

where Vj = j-dj, p^ = p^ + 5a39i(/(a;j))/(2|/p) and T"^'^ is a completely antisymmetric 
tensor with -T^^s ^ 2^231 ^ 2"3i2 = 1, 

The solution to this equation is as follows 

p': = R^'^'q^V^q, (3.9) 

where qp satisfy qi + q2 ~ q^ ~ qi = 1 and i?"^"" is a tensor, antisymmetric in the last 
two indices, with i^^^i ^ ^i34 ^ ^213 ^ ^242 ^ ^3i4 ^ ^323 ^ ^^ rj^^^ functions q^ then 

establish a mapping from S'f x K onto that hyperboloid. 

If we now plug the solution for p' Eq.(3.9) into the self-duality equation, and after 
considerable massaging, we arrive at 

5os" = e"'^Sb4 (3.10) 

Vis', = VjVkSa {^i^J^k) (3.11) 

s'^sa = l/T (3.12) 

where indices are raised and lowered with the metric gab = diag{l^— 1,-1) and Si = 
1/1(9? + <?! + ?| + qI), S2 = 2|/|(g2g3 - qiq4) and S3 = 2|/|(gig3 + g2g4). The dynamics 
Eq.(3.10) is consistent with the constraint Eq.(3.12). The most important restriction 
follows from the integrability condition of Eq.(3.11), which together with Eq.(3.12) serves 
to restrict the possible value of f{x). If we introduce the function Zi = z{xi) by the 
condition -^^ = f{x), Eq.(3.11) forces Sa to be sum of single variable functions of the 
variables 01^1 + 02-22 + 03^3 with Oj = +1, —1 or and J^o-i = ^{iTLod2). On the other hand 
l/P is a sum F{zi) + F{z2) + -^(-23)- One class of solutions is given by functions Sa{zi + 
-22 + -23)- In this case Eq.(3.12) forces f{x) = Ax and then the whole problem possesses 
spherical symmetry Zi + Z2 + z^ oc |xp. We are then in the situation discused in Ref.9 
and our equations recover all the axially symmetric solutions including the B.P.S.T.[10] 
one. It is not easy to find other solutions of Eq.(3.11) and Eq.(3.12). Indeed, if we impose 
these conditions on the coefficients of the Taylor expansion of the functions Sa{z), the 
number of conditions grows faster than the number of parameters. This suggests to try 
polynomials for the functions Sa- In this way we have discovered other solutions, most 
remarkably one which is defined on the torus with z = Acos{Bx + C) and where si and 
S3 are polynomials of second degree in Zi and S2 of first degree. The solution goes to a 
pure gauge in t ^ ±00 but is singular at t = and the total action is divergent ( in the 
unit cell ). We have been unable to find any solution which is satisfactory and coincides 
with our numerical data. 
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If our solution does not follow from Eqs.(3.10)-(3.12), one has to drop some of the 
assumptions leading to those equations. Most likely it is Eq.(3.6) with the assumed 
properties for p^ which is wrong, since its validity is simply based on the ability to describe 
the data. Unfortunately, if we go back to the form Eq.(3.1) with Q*-^-* = Q^^^-* + Q^^~^ 
we have no additional equations which could allow an analytic development of the sort 
leading to Eqs.(3.10)-(3.12). One can simply rewrite the self duality equations for the Q 
functions, but we do not know how to select solutions going to a pure gauge in t = ±00. 

4 Summary and conclusions 

In this paper, we have studied the instanton-like solution which occurs in a 3-torus with 
symmetric twist. It is remarkable that the cooling method provides a very accurate 
description of this solution which allows a systematic analysis of its properties. Our re- 
sults imply that there are indeed four instanton and four anti-instanton solutions, all 
of which are (within errors) cubic symmetric and P ■ T invariant. The solution is also 
smooth, although singularities in sufficiently high derivatives are not excluded. We have 
given functional forms which incorporate the previously mentioned properties, plus ad- 
ditional restrictions consistent with the data. Nonetheless, the more restrictive form, 
which allows the reduction of the problem to motion on a hyperboloid and recovers many 
multi-instanton solutions on S4 does not seem to contain the solution of our problem. 
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Table Captions 

Table 1: The value of the parameters obtained when fitting the functional form (2.11) and 
(3.2) to the Ng = 11, Nf = 29 data for A\. rit = 15 corresponds to the time of 
maximum energy. Errors are indicated by giving within parenthesis their magnitude 
affecting the last quoted digit. 

Table 2: The same as Table 1 but for the parameters entering in Af, A^. 

Table 3: The same as Table 1 but for the data of -Bj. 

Table 4: The same as Table 2 but for the data of Bf, Bf. 
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Figure Captions 

Fig.l: The scalar product M12 = Ei ■ E2 plotted as a function of x for y = 0.209 and 
z = 0. The squares correspond to the lattice values for iV^ = 11, A^j = 21. They are 
compared to the prediction of Eq.(2.3). 

Fig.2: rrioit) appearing in Eq(2.5) is plotted as a function of time for Ng = 11, Nt = 21. 
The continous line represents the function cos(|tan/i(0.697rt)). 
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Table 1 



nt 


qi 


q2 


qs 


q4 


5 


0.094(8) 


0.04(6) 


0.02(6) 


0.00(1) 


6 


0.140(8) 


0.05(2) 


0.02(2) 


0.00(1) 


7 


0.208(8) 


0.06(3) 


0.02(3) 


0.00(1) 


8 


0.302(8) 


0.07(2) 


0.02(2) 


0.00(1) 


9 


0.438(8) 


0.09(1) 


0.03(2) 


0.00(1) 


10 


0.624(8) 


0.103(9) 


0.043(9) 


0.00(1) 


11 


0.870(8) 


0.120(6) 


0.058(7) 


0.01(1) 


12 


1.182(6) 


0.138(4) 


0.078(5) 


0.02(1) 


13 


1.550(6) 


0.154(3) 


0.102(4) 


0.02(1) 


14 


1.954(6) 


0.167(3) 


0.129(3) 


0.03(1) 


15 


2.356(6) 


0.177(1) 


0.156(2) 


0.05(1) 


16 


2.716(6) 


0.182(2) 


0.182(2) 


0.07(1) 


17 


3.010(6) 


0.184(1) 


0.205(2) 


0.09(1) 


18 


3.234(6) 


0.183(1) 


0.224(2) 


0.11(1) 


19 


3.392(6) 


0.181(1) 


0.240(2) 


0.12(1) 


20 


3.498(6) 


0.179(1) 


0.252(2) 


0.13(1) 


21 


3.566(6) 


0.178(1) 


0.261(2) 


0.14(1) 


22 


3.610(6) 


0.176(1) 


0.267(2) 


0.14(1) 


23 


3.638(6) 


0.175(1) 


0.272(2) 


0.14(1) 


24 


3.654(6) 


0.174(1) 


0.276(2) 


0.14(1) 


25 


3.664(6) 


0.174(1) 


0.278(2) 


0.14(1) 


26 


3.672(6) 


0.173(1) 


0.280(2) 


0.14(1) 


27 


3.674(6) 


0.173(1) 


0.282(2) 


0.14(1) 


29 


3.674(6) 


0.172(1) 


0.284(2) 


0.14(1) 
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Table 2 



nt 


ai 


a2 


as 


a4 


as 


6 


0.00(6) 


-0.02(6) 


-0.122(6) 


0.07(4) 


-0.002(6) 


7 


0.01(5) 


-0.03(4) 


-0.172(6) 


0.09(3) 


0.008(6) 


8 


0.01(3) 


0.01(3) 


-0.244(6) 


0.11(2) 


0.018(6) 


9 


0.01(2) 


0.03(2) 


-0.348(6) 


0.14(2) 


0.034(6) 


10 


0.02(1) 


0.04(1) 


-0.496(6) 


0.15(1) 


0.062(6) 


11 


0.02(8) 


0.07(1) 


-0.704(6) 


0.19(8) 


0.109(6) 


12 


0.033(6) 


0.092(5) 


-0.984(6) 


0.225(6) 


0.174(6) 


13 


0.042(5) 


0.118(5) 


-1.344(6) 


0.253(4) 


0.272(6) 


14 


0.052(4) 


0.145(4) 


-1.776(6) 


0.278(3) 


0.402(6) 


15 


0.060(3) 


0.169(3) 


-2.258(6) 


0.298(3) 


0.558(6) 


16 


0.065(3) 


0.189(2) 


-2.746(6) 


0.311(2) 


0.730(6) 


17 


0.068(2) 


0.204(2) 


-3.200(6) 


0.314(2) 


0.902(6) 


18 


0.069(2) 


0.215(2) 


-3.590(6) 


0.324(2) 


1.056(6) 


19 


0.069(2) 


0.221(2) 


-3.902(6) 


0.325(2) 


1.186(6) 


20 


0.067(2) 


0.225(1) 


-4.138(6) 


0.325(1) 


1.290(6) 


21 


0.065(2) 


0.227(2) 


-4.312(6) 


0.324(1) 


1.368(6) 


22 


0.064(2) 


0.229(1) 


-4.434(6) 


0.323(1) 


1.424(6) 


23 


0.063(2) 


0.230(1) 


-4.520(6) 


0.322(1) 


1.436(6) 


24 


0.062(2) 


0.230(1) 


-4.580(6) 


0.321(1) 


1.494(6) 


25 


0.061(2) 


0.231(1) 


-4.620(6) 


0.320(1) 


1.514(6) 


26 


0.061(2) 


0.231(1) 


-4.648(6) 


0.320(1) 


1.530(6) 


27 


0.061(2) 


0.231(1) 


-4.666(6) 


0.320(1) 


1.540(6) 


29 


0.061(2) 


0.232(1) 


-4.686(6) 


0.320(1) 


1.554(6) 
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Table 3 



nt 


qi 


q2 


qs 


q4 


3 


0.216(8) 


0.04(3) 


0.00(3) 


0.00(1) 


5 


0.424(8) 


0.06(1) 


0.01(2) 


0.00(1) 


6 


0.610(8) 


0.07(8) 


0.02(1) 


0.00(1) 


7 


0.876(8) 


0.090(5) 


0.026(8) 


0.00(1) 


8 


1.246(8) 


0.109(4) 


0.038(6) 


0.01(1) 


9 


1.736(8) 


0.130(3) 


0.054(4) 


0.01(1) 


10 


2.352(8) 


0.152(2) 


0.077(3) 


0.03(1) 


11 


3.058(8) 


0.175(2) 


0.107(2) 


0.05(1) 


12 


3.760(8) 


0.196(1) 


0.144(2) 


0.08(1) 


13 


4.300(8) 


0.213(1) 


0.191(2) 


0.12(1) 


14 


4.504(8) 


0.223(1) 


0.244(2) 


0.17(1) 


15 


4.278(6) 


0.224(1) 


0.301(2) 


0.21(1) 


16 


3.682(8) 


0.216(1) 


0.361(2) 


0.23(1) 


17 


2.900(8) 


0.200(2) 


0.421(2) 


0.21(1) 


18 


2.124(8) 


0.181(2) 


0.482(3) 


0.17(1) 


19 


1.474(6) 


0.160(3) 


0.542(5) 


0.12(1) 


20 


0.986(6) 


0.141(4) 


0.603(7) 


0.08(1) 


21 


0.644(6) 


0.124(6) 


0.663(1) 


0.04(1) 


22 


0.412(6) 


0.108(8) 


0.73(1) 


0.03(1) 


23 


0.264(6) 


0.09(1) 


0.79(3) 


0.01(1) 


24 


0.166(6) 


0.08(2) 


0.87(4) 


0.00(1) 


26 


0.062(6) 


0.03(4) 


1.1(1) 


0.00(1) 
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Table 4 



nt 


ai 


a2 


as 


a4 


as 


3 


0.00(6) 


0.00(6) 


-0.120(6) 


0.07(4) 


0.006(6) 


5 


0.01(2) 


0.02(3) 


-0.292(6) 


0.10(2) 


0.022(6) 


6 


0.01(2) 


0.02(2) 


-0.438(6) 


0.12(1) 


0.042(6) 


7 


0.02(1) 


0.04(1) 


-0.644(4) 


0.147(8) 


0.074(6) 


8 


0.019(2) 


0.050(8) 


-0.940(6) 


0.177(6) 


0.132(6) 


9 


0.027(5) 


0.069(5) 


-1.354(6) 


0.210(4) 


0.224(6) 


10 


0.037(4) 


0.095(4) 


-1.912(3) 


0.245(3) 


0.370(6) 


11 


0.050(3) 


0.126(3) 


-2.626(6) 


0.280(2) 


0.584(6) 


12 


0.064(2) 


0.161(2) 


-3.464(6) 


0.313(2) 


0.872(6) 


13 


0.079(2) 


0.198(1) 


-4.322(6) 


0.341(1) 


1.212(6) 


14 


0.092(2) 


0.233(1) 


-5.018(6) 


0.361(1) 


1.550(6) 


15 


0.100(2) 


0.262(1) 


-5.354(6) 


0.372(1) 


1.796(6) 


16 


0.100(2) 


0.282(1) 


-5.214(6) 


0.372(1) 


1.880(6) 


17 


0.093(2) 


0.293(1) 


-4.754(6) 


0.364(1) 


1.786(6) 


18 


0.080(2) 


0.297(2) 


-3.950(6) 


0.349(2) 


1.556(6) 


19 


0.064(3) 


0.294(2) 


-2.998(6) 


0.331(2) 


1.264(6) 


20 


0.047(4) 


0.289(3) 


-2.232(6) 


0.312(3) 


0.974(6) 


21 


0.031(5) 


0.282(4) 


-1.608(6) 


0.293(4) 


0.724(6) 


22 


0.017(7) 


0.276(6) 


-1.132(6) 


0.277(3) 


0.522(6) 


23 


0.00(1) 


0.271(8) 


-0.784(6) 


0.264(8) 


0.370(6) 


24 


0.00(1) 


0.27(1) 


-0.536(3) 


0.26(1) 


0.230(6) 


25 


0.00(2) 


0.28(1) 


-0.364(6) 


0.26(2) 


0.180(6) 


26 


0.00(3) 


0.30(2) 


-0.244(6) 


0.27(3) 


0.130(6) 


29 


0.00(3) 


0.30(2) 


-0.244(6) 


0.27(3) 


0.130(6) 
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